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INTRODUCTION 


The  primary  objective  of  this  research  centered  around  breast  cancer  imaging,  specifically 
developing  new  techniques  in  single  photon  emission  computed  tomography  (SPECT)  breast  cancer 
imaging  which  will  provide  improved  detection  and  characterization  of  breast  cancer.  The  leading 
hypothesis  in  this  research  is  that  by  using  a  vertical  axis  of  rotation  (VAOR)  and  various  camera 
head  and  collimator  configurations,  several  advantages  can  be  gained  over  traditional  horizontal 
axis  of  rotation  (HAOR)  SPECT.  A  tiltable  head  SPECT  (TH-SPECT)  system  allows  for  improved 
imaging  of  lesions  closer  to  the  chest  wall  as  well  as  a  smaller  radius  of  rotation  [1].  Compared  to  a 
filtered  backprojection  reconstruction  algorithm,  OSEM  reconstruction  algorithm  displayed 
improved  images  [3].  The  use  of  a  pinhole  collimator  allows  for  less  background  activity  to  be 
detected,  a  smaller  radius  of  rotation,  and  higher  values  for  contrast  and  SNR  [4]. 


BODY 


Disadvantages  of  HAOR  breast  imaging  with  SPECT  include  attenuation  through  the  torso 
and  a  large  radius  of  rotation  (ROR).  Tilting  a  camera  head  allows  the  collimator  to  get  closer  to  the 
breast  and  chest  wall.  Using  various  tilt  angles  for  the  camera  heads,  an  axially  Tiltable  2-Headed 
SPECT  (TH-SPECT)  system  was  used  to  observe  resolution,  contrast  and  SNR  values  in  tillable 
breast  and  torso  phantoms  containing  a  1  cm  lesion  suspended  in  the  breast.  The  results  of  this 
HAOR  TH-SPECT  system  were  compared  to  high-count  planar  images.  An  ordered  subsets 
expectation  maximization  (OSEM)  algorithm  (8  subsets,  3  iterations),  which  accounted  for  the 
geometry  of  the  system,  was  used  to  reconstruct  the  data.  In  comparison  to  the  high  count  planar 
images,  a  0.5mL  lesion  was  more  clearly  visible  in  TH-SPECT.  Lesion  contrast  and  lesion  SNR 
improved  nearly  three  and  two  times,  respectively.  As  a  result,  this  system  configuration  reduced 
the  radius  of  rotation  and  allowed  for  adequate  imaging  near  the  chest  wall  breast  tissue,  improving 
sensitivity  and  resolution  [  1  ] .  A  more  detailed  description  of  these  results  can  be  found  in  the 
appendix  in  “Breast  Tumor  Imaging  Using  Tiltable  Head  SPECT  Cameras.”  The  one  drawback  to 
this  system  is  the  fact  that  axially  blurring  affects  the  chest  wall  region  of  the  breast.  Future  work 
could  be  done  with  a  combination  of  two  heads  tilting  at  different  angles,  aiming  to  provide  more 
complete  sampling. 

One  concern  dealing  with  TH-SPECT  that  was  evident  in  the  previous  study  dealt  with 
insufficient  axial  sampling  of  the  image  space.  Specifically,  insufficient  sampling  of  the 
myocardium,  liver  and  torso  will  cause  interference  when  imaging  small  lesions  in  the  breast  and 
thus  result  in  artifacts  in  the  reconstructed  images.  Tilting  the  camera  head  prevents  the 
reconstruction  algorithms  from  meeting  the  Orlov  sphere  requirement,  which  states  that  an  image 
space  must  be  sampled  over  a  great  circle  [2].  A  three-dimensional  algorithm  is  needed  to 
reconstruct  TH-SPECT  data.  Iterative  methods,  such  as  OSEM,  have  been  used  to  model  the 
sampling  of  frequency  space  associated  with  circular  orbits  with  a  tilted  geometry.  In  the  following 
research,  a  filtered  backprojection  (FBP)  algorithm  was  devised  with  the  intent  of  successfully 
reconstructing  tilted-head  geometries.  This  FBP  method  was  compared  to  OSEM  for  a  fillable 
breast  phantom  containing  lesions  at  various  tilt  angles  and  for  a  fillable  breast  phantom  containing 
a  mini-Defrise  disk  imaged  with  parallel-beam  TH-SPECT  at  angles  of  0-15  degrees.  Values  for 
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both  contrast  and  SNR  were  compared.  The  publication  supporting  this  research  can  be  found  in 
the  appendix.  In  this  paper  [3],  a  geometric  derivation  of  the  ramp  filter  is  presented  for  tilted 
parallel  beam  geometries.  Using  this  ramp  filter,  a  FBP  method  was  implemented  and  compared 
with  an  iterative  OSEM  algorithm  using  TH-SPECT  data.  The  FBP  results  indicated  a  greater 
increase  in  backprojection  artifacts  with  greater  angular  tilt,  as  well  as  a  greater  bias  in  FBP  noise 
versus  bias  plots  with  a  greater  angular  tilt.  OSEM  reconstructions  showed  greater  SNR  and 
contrast  values  at  all  tilt  angles  [3].  Because  TH-SPECT  data  do  not  sufficiently  sample  the  Orlov 
sphere,  artifacts  are  present  in  all  reconstructions.  However,  OSEM  TH-SPECT  reconstruction 
offers  better  shape  and  uniform  activity  distribution  of  the  breast  compared  to  FBP  methods  [3]. 

A  pinhole  collimator  is  advantageous  over  a  parallel-hole  collimator  in  imaging  the  breast 
in  its  ability  to  decrease  the  ROR,  resulting  in  better  resolution  for  detecting  lesions.  However, 
an  issue  in  both  parallel  and  pinhole  SPECT  breast  imaging  is  photon  contamination  from  other 
organs  that  will  impact  the  clarity  of  the  image,  such  as  the  liver,  myocardium,  lungs  and  gall 
bladder.  Since  a  pinhole  collimator,  due  to  its  geometry,  will  collect  fewer  photons  from  the 
torso,  the  hypothesis  that  we  set  out  to  prove  was  that  lead  shielding  on  the  bed  would  further 
reduce  the  unwanted  photons  from  outlying  organs  that  were  collected  by  the  pinhole  collimator 
and  thus  increase  the  contrast  and  SNR  of  the  image.  Using  a  Trionix  XLT  scanner  (Trionix 
Research  Lab,  Twinsburg,  Ohio)  with  one  pinhole  collimator,  experimental  data  was  obtained  of 

a  breast  phantom  mounted  on  an  anthropomorphic 
torso  phantom  placed  in  a  prone  position  on  the  bed 
(Figure  1).  The  first  set  of  scans  consisted  of  the 
torso  phantom  and  all  its  organs  filled  with  activity, 
and  the  second  set  of  scans  consisted  of  the  same 
configuration  of  the  torso  phantom  as  the  first  scan, 
except  with  the  addition  of  lead  shielding  placed  over 
the  entire  table  except  for  where  the  breast  was  located.  The  angular  range  of  these  scans  was 
186  degrees,  resulting  in  an  incomplete  circular  orbit  (ICO).  The  breast  contained  two  spherical 
lesions  (1.6  cm  diameter),  both  of  which  were  filled  with  activity.  The  projection  data  were 
reconstructed,  without  shielding  attenuation  modeling  or  scattering  and  attenuation  correction, 
with  an  OSEM  algorithm  having  thirty  iterations  and  eight  subsets.  A  3-D  region-of-interest 
analysis  was  performed  on  the  lesions  to  evaluate  shielding  effects.  The  use  of  shielding  slightly 
reduced  noise  in  contrast  estimates,  from  about  0.25  to  0.23.  However,  shielding  increased  bias 
in  contrast  from  0.24±0.08  to  0.44±0.07.  SNRs  were  also  lower  with  shielding.  These  results 
suggest  that  if  shielding  is  to  be  of  benefit,  shielding  attenuation  must  be  modeled  within  iterative 
reconstruction.  However,  the  standard  reconstruction  algorithms  such  as  OSEM  and  MLEM  are 
limited  in  reconstructing  inconsistent  projection  data.  An  abstract  [4]  and  supplementary 
information  was  submitted  accepted  to  the  2004  NSS/MIC  conference  and  is  also  attached  for  a 
more  detailed  explanation  of  this  experiment. 

■  An  advantage  of  the  Summit  T-22  SPECT  system,  shown  in  Figure  2,  is 

the  fact  that  the  camera  heads  can  be  tilted  through  varying  angles.  Currently, 
the  T-22  can  only  use  parallel  collimators.  To  tie  together  all  the  previous 
research,  we  want  to  use  a  tiltable  head  system  with  a  small  radius  of  rotation 
in  conjunction  with  a  pinhole  collimator  to  image  the  breast.  In  order  to  use 
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Figure  2 
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the  pinhole  collimator  in  conjunction  with  the  T-22  camera,  we  have  designed  an  adaptor 
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Figure  3 


plate,  shown  in  Figure  3,  to  allow  the  pinhole  collimator  to 
be  mounted  onto  the  camera.  This  plate  was  designed 
using  Computer  Assisted  Design  and  is  in  the  process  of 
being  built  by  the  Duke  University  Physics  Shop.  There 
are  actually  two  plates  being  built,  one  for  each  of  the 
detector  heads  on  the  T-22  system.  Once  these  plates  are 
finished,  the  pinhole  collimator  can  be  attached  to  the  T- 
22  system  and  used  to  image  the  breast.  For  the 
acquisition  of  experimental  VAOR  phantom  studies,  two 
pinhole  collimators  would  be  mounted  on  the  two  tilted 
detector  heads.  A  different  tilt-angle  would  be  used  for 
each  detector  head. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Development  and  evaluation  of  FBP  TH-SPECT  image  reconstruction  algorithm  using  a 
geometric  derivation  of  the  ramp  filter  [3]. 

•  Comparison  of  FBP  TH-SPECT  algorithm  with  traditional  FBP  and  OSEM  [3]. 

•  Acquisition,  evaluation,  and  comparison  of  TH-SPECT  data,  including  contrast  and  S/N 
analysis  using  experimental  phantoms  [1], 

•  Acquisition  and  evaluation  of  combined  horizontal  and  vertical  axes  of  rotation. 

•  Comparison  of  contrast  and  SNR  values  between  shielded  and  unshielded  torso  phantom  [4]. 

•  Design  of  special  adapter  plate  for  using  Triad  XLT  pinhole  collimator  on  T22  SPECT 
system. 


REPORTABLE  OUTCOMES: 

Manuscripts 

•  Pieper  BC,  Bowsher  JE,  Tomai  MP,  Archer  CN,  Jaszczak  JJ,  Parallel-Beam  Tilted-Head  Analytic 
SPECT  Reconstruction:  Derivation  and  Comparison  with  OSEM,  IEEE  Trans.  Nucl.  Sci.,  vol.  49, 
pp.  2394-2400,  Oct  2002. 

•  Pieper  BC,  Bowsher  JE,  Tomai  MP,  Peter  J,  Jaszczak  JJ,  Breast  Tumor  Imaging  Using  Tiltable  Head 
SPECT  Cameras, 

•  Lee  T,  Braun  K,  Jazczak  J,  Bowsher  J,  Bobkov  K,  Effects  of  Lead  Shielding  of  Background 
Organ  Activity  in  Pinhole  SPECT  Breast  Tumor  Imaging,  Submitted  to  the  2004  IEEE 
Nuclear  Science  Symposium  &  Medical  Imaging  Conference,  16-21  October,  2004,  Rome, 
Italy. 
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CONCLUSIONS 


>> 


From  the  above  research,  several  conclusions  can  be  drawn.  To  combat  the  issue  of 
insufficient  sampling,  it  was  shown  that  an  OSEM  algorithm  will  reconstruct  the  projection  data 
better  than  a  FBP,  showing  greater  SNR  and  contrast  values  at  all  tilt  angles.  To  combat  the 
issues  produced  by  a  HAOR,  the  camera  heads  were  tilted,  reducing  the  radius  of  rotation  and 
allowing  for  adequate  imaging  near  the  chest  wall  breast  tissue.  The  combination  of  an  OSEM 
algorithm  and  TH-SPECT  data  offers  better  shape  and  uniform  activity  distribution  of  the  breast, 
improving  contrast  and  lesion  SNR.  Concerning  the  lead  shielding,  comparison  of  contrast  and 
SNR  between  shielded  and  non-shielded  images  revealed  that  contrast  and  SNR  is  significantly 
lower  with  lead  shielding  than  without  lead  shielding,  the  opposite  of  what  we  had  expected. 

In  the  future,  applications  of  SPECT  to  breast  imaging  may  provide  an  alternative  method 
of  detecting  breast  cancer.  Future  research  is  required  to  determine  whether  an  OSEM  algorithm 
can  be  implemented  that  can  effectively  account  for  inconsistent  projection  data.  Also,  future 
work  may  include  developing  more  elaborate  orbits  that  provide  a  complete  sampling  of  the 
breast,  thus  meeting  the  Orlov  criteria. 
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Breast  Tumor  Imaging  Using  a  Tillable  Head  SPECT 

Camera 
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Jorg  Peter,  Member,  IEEE ,  Kim  Greer,  and  Ronald  J.  Jaszczak,  Fellow,  IEEE 


Abstract — As  a  test  bed  for  dedicated  breast  single-photon 
emission  computed  tomography  (SPECT)  cameras  that  are  under 
development,  a  general  purpose  SPECT  system  with  two  tiltable 
heads  was  used  to  image  tillable  breast  and  torso  phantoms  con¬ 
taining  multiple  lesions.  Breast,  liver,  and  myocardial  activity  were 
included  in  order  to  simulate  direct  contamination  and  Compton 
scattering  expected  in  clinical  scans.  The  tiltable-head  SPECT 
(TH-SPECT)  data  were  reconstructed  using  an  OS-EM  algorithm 
which  accounted  for  the  tilted  geometry.  High-count  planar 
images  were  acquired  for  comparison  with  TH-SPECT.  In  order 
to  characterize  axial  blurring  effects  inherent  with  TH-SPECT 
reconstructions,  two  cylindrical  disk  Defrise  phantoms,  one  large 
Defrise  phantom  and  one  mini-Defrise  phantom  placed  inside 
the  fillable  breast  phantom,  were  imaged  at  various  tilt  angles. 
Results  indicate  an  increase  in  axial  blurring  with  greater  tilt 
angle.  Reconstructions  of  the  combined  fillable  breast  and  torso 
phantoms  containing  two  1.15  ml  lesions,  one  centered  axially 
and  one  proximal  to  the  anterior  chest  wall  within  the  breast, 
were  most  clearly  visible  in  the  30°  reconstructed  TH-SPECT 
images,  providing  lesion  signal-to-noise  ratio  (SNR)  and  contrast 
improvements  of  nearly  three  times  compared  to  the  high-count 
planar  images. 

I.  Introduction 

FOR  BREAST  cancer  imaging  using  conventional  single¬ 
photon  emission  computed  tomography  (SPECT)  systems, 
the  gamma  camera’s  axis  of  rotation  is  parallel  to  the  patient 
table— a  horizontal  axis  of  rotation  (HAOR).  This  geometry  de¬ 
creases  breast  signals  due  to  torso  attenuation  and,  more  im¬ 
portantly,  requires  a  large  radius  of  rotation  which  yields  poor 
results  for  SPECT  breast  imaging  [1]. 

Previous  work  indicates  breast  cancer  imaging  can  be 
improved  by  using  a  vertical  axis  of  rotation  (VAOR)  [1],  [2]. 
In  this  paper,  we  expand  upon  previous  VAOR  breast  imaging 
by  using  tiltable  gamma  cameras.  By  tilting  the  camera  off 
axis  with  respect  to  the  nipple-to- anterior  chest  axis  through 
the  breast,  improved  imaging  of  lesions  near  the  chest  wall 
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Fig.  1.  Schematic  of  dedicated  breast  VAOR  system. 

and  a  reduced  radius  of  rotation  (ROR)  can  be  achieved.  This 
paper  investigates  tilted  VAOR  breast  imaging  (see  Fig.  1)  by 
simulating  a  dedicated  VAOR  system  using  a  HAOR  SPECT 
system  with  tiltable  gamma  cameras. 

One  major  concern  with  tiltable-head  SPECT  (TH-SPECT) 
imaging  is  insufficient  angular  sampling.  The  completeness  of 
reconstructions  using  TH-SPECT  is  bound  by  Orlov  sphere  con¬ 
ditions  [3].  In  order  to  satisfy  these  conditions,  the  image  space 
must  be  sampled  over  a  great  circle  or  equatorial  orbit.  By  tilting 
the  camera  heads  off  axis,  this  primary  requirement  is  not  met 
(see  Fig.  2),  and  thus  there  is  the  potential  for  artifacts  in  the 
reconstructed  images.  This  is  also  an  issue,  in  general,  with 
slant-hole  collimators  and  rotating  slant-hole  collimators  [4]. 

Another  concern  with  TH-SPECT  imaging  is  background  in¬ 
terference  from  myocardium,  liver,  and  body  tissue  that  is  ex¬ 
pected  clinically  with  the  use  of  Tc-99m  labeled  tracers.  By 
tilting  the  camera  heads  off  axis,  heart  and  liver  activity  are  in¬ 
troduced  into  the  camera’s  field  of  view  (FOV)  at  some  projec¬ 
tions,  potentially  contaminating  TH-SPECT  reconstructions. 

II.  Materials  and  Methods 

All  data  were  acquired  using  a  tiltable  two-headed  SPECT 
system  (Vision  T-22,  Summit  Nuclear,  Twinsburg,  OH)  (see 
Fig.  5).  Only  one  head  was  used  to  obtain  data  during  this  study. 
A  parallel  beam  low-energy  high-resolution  (LEHR)  collimator 
(Hitachi)  was  used  in  all  experiments.  Collimator  hole  length 
was  41.0  mm,  hexagonal  hole  size  was  1.8  mm,  flat-to-flat  and 
septal  thickness  was  0. 1 8  mm.  All  images  were  reconstructed  by 
OS-EM  (5  iterations,  8  subsets),  with  the  tilted  geometry  mod¬ 
eled  in  the  photon  detection  probabilities.  Compton  scattering, 
attenuation,  and  collimator  and  crystal  spatial  resolution  were 
not  modeled. 
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Fig.  4.  Photograph  of  the  custom  breast  phantom  with  three  mini-Defrise  disks 
8.0  cm  in  diameter  and  1 .0  cm  thick,  separated  by  a  1 .0  cm  gap  and  aligned  along 
the  breast  axis.  Nipple  is  at  right,  and  curvature  at  left  corresponds  with  torso 
shape. 


Fig.  3.  Schematic  of  the  Defrise  phantom  with  camera  head  at  0°,  15° ,  30°, 
45° ,  and  60°  tilt  angles.  Curved  arrow  illustrates  camera  orbit  out  of  page  about 
the  AOR. 

A.  Defrise  Disk  Phantom 

A  Defrise  disk  phantom  (see  Fig.  3)  (Data  Spectrum  Corp., 
Hillsborough,  NC)  was  utilized  in  the  first  experiment.  This 
phantom  consisted  of  a  cylinder  containing  six  acrylic  disks 
placed  1.5  cm  apart.  The  acrylic  disks  had  a  thickness  of  1.2 
cm  and  a  diameter  of  18.9  cm.  Five  emission  disks  were  then 
formed  by  adding  1 6.6  mCi  of  Tc-99m  and  water  to  the  cylinder, 
each  hot  disk  having  a  uniform  concentration  of  4.1  /rCi/ml. 
The  phantom  was  then  placed  parallel  to  the  camera’s  axis  of 
rotation  (AOR)  with  a  radius  of  rotation  (ROR)  of  11.5  cm, 
measured  from  the  center  of  the  camera  to  the  Defirse  disk 
phantom’s  long  axis,  perpendicular  to  the  camera’s  AOR  (see 
Fig.  3).  Projection  data  were  taken  for0°,  15°,  30°,  45°,  and  60° 
tilt  angles.  Each  projection  data  set  was  acquired  for  128  angular 
views  over  360°,  and  each  projection  was  128  x  128  pixels  and 
pixel  size  was  0.48  cm  square.  Total  acquisition  time  was  7  s  per 
angular  step  for  a  total  time  of  15  minutes.  Subsequent  acquisi¬ 
tion  times  were  increased  to  compensate  for  radioactive  decay. 

B.  Breast  Phantom  with  Mini-Disk  Phantom 

In  order  to  characterize  axial  blurring  within  the  breast  geom¬ 
etry,  custom  mini-Defrise  phantom  disks  (Data  Spectrum  Corp., 
Hillsborough,  NC)  were  placed  inside  a  custom  fillable  breast 
phantom  (Data  Spectrum  Corp.,  Hillsborough,  NC)  (see  Fig.  4). 
By  filling  the  breast  phantom  with  water  having  a  22.5  //.Ci/mL 


concentration  of  Tc-99m,  three  cold  disks  measuring  1.0  cm 
in  thickness  and  8.0  cm  in  diameter  and  two  hot  disks  mea¬ 
suring  1.0  cm  and  8.0  cm  in  diameter  were  formed.  The  axis 
of  the  breast  phantom  was  then  placed  parallel  to  the  camera’s 
AOR  and  the  radius  of  rotation  was  measured  from  the  nipple 
to  the  center  of  the  camera,  perpendicular  to  the  camera’s  AOR. 
Because  axial  blurring  becomes  extremely  noticeable  at  larger 
tilt  angles  (>  15°),  projection  data  were  taken  for  angles  of 
0°,  5°,  10°,  and  15°  having  RORs  of  12.1  cm,  11.8  cm,  11.6 
cm,  and  10.9  cm  respectively.  The  ROR  was  decreased  with 
greater  tilt  angle  due  to  greater  access  to  the  breast.  Potentially 
smaller  RORs  could  be  possible  with  dedicated  small  area  breast 
imaging  systems  [5].  Total  acquisition  time  was  10  s  per  angular 
step  or  21.33  min.  Subsequent  acquisition  times  were  increased 
to  compensate  for  radioactive  decay. 

C.  Anthropomorphic  Torso  Phantom  and  Breast  Phantom 
with  Lesions 

An  anthropomorphic  torso  phantom  (Model  ECT/TOR/P, 
Data  Spectrum  Corp.,  Hillsborough,  NC)  combined  with  the 
fillable  breast  phantom  was  imaged  in  this  experiment  (see 
Fig.  5)  to  simulate  the  expected  emission  contamination  from 
the  myocardium  and  liver  to  the  specific  breast  imaging  task. 
The  craniocaudal  axis  of  the  torso  phantom  was  perpendicular 
to  the  camera’s  axis  of  rotation,  thus  necessitating  the  tilted 
camera  heads  and  allowing  for  close  proximity  imaging  of  the 
breast  phantom.  Two  1.3  cm  diameter  lesions  (1.15  ml)  with 
Tc-99m  concentration  of  11.1  /xCi/ml  were  suspended  inside 
the  left  breast.  One  lesion  was  placed  axially  in  the  center  of  the 
breast,  5.1  cm  from  the  nipple..  The  second  lesion  was  located 
2.2  cm  anterior  to  the  chest  wall  and  positioned  3.0  cm  toward 
the  medial-superior  side  of  the  breast.  A  third  0.9  cm  diameter 
lesion  (0.38  ml)  with  a  concentration  of  11.1  //Ci/ml  was 
placed  outside  the  breast  on  the  chest  wall  to  represent  axillary 
tumor  growth.  The  breast  phantom  had  a  uniform  concentration 
of  1.5  /xCi/ml.  In  order  simulate  contamination  effects  due  to 
direct  superposition  of  myocardial  activity  on  the  breast  at 
some  projection  views  and  also  Compton  scattering  expected 
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Fig.  5.  Photograph  of  anthropomorphic  torso  and  breast  phantoms  with  dual 
tilted  heads  on  the  SPECT  camera. 


in  clinical  scans,  the  tumor  and  breast  were  imaged  with 
additional  liver  and  myocardium  activity.  The  lesion  diver:  my¬ 
ocardium  rbreast-background  ratio  was  7.4:12.5:12.5:1,  giving 
the  liver  and  heart  a  Tc-99m  concentration  of  18.7  //Ci/ml. 

Projection  data  were  acquired  for  0°,  15°,  30°,  and  45°  tilt 
angles  with  ROR’s  of  18.1  cm,  10.5  cm,  9.4  cm,  and  5.9  cm  re¬ 
spectively,  measured  as  previously  described.  The  images  were 
acquired  for  128  angular  views  over  360°,  and  each  projection 
was  128  x  128  pixels  and  pixel  size  was  0.48-cm  square.  Total 
acquisition  time  was  32  min,  or  15  s  per  angular  step  and  sub¬ 
sequent  acquisition  times  were  increased  to  compensate  for  ra¬ 
dioactive  decay. 

Because  planar  scintimammography  is  gaining  acceptance  as 
an  effective  way  to  identify  breast  lesions  greater  than  1  cm  in 
diameter  [6],  a  planar  scintimammography  image  (craniocaudal 
view)  was  made  for  comparison.  Breast,  liver,  and  myocardium 
background  radiation  were  present,  with  the  same  activity  con¬ 
centration  as  for  TH-SPECT,  and  the  total  acquisition  time  was 
5  min. 


III.  Results 


Fig.  6.  Single  slice  sagittal  views  of  OS -EM  reconstructed  Defrise  phantom 
images  at  indicated  tilt  angles. 


A.  Degradative  Angular  Sampling  Effects 

Single  slice  sagittal  views  of  the  reconstructed  Defrise 
phantom  images  at  0°,  15°,  30°,  45°,  and  60°  are  presented  in 
Fig.  6.  For  the  reconstructed  sagittal  slices,  each  disk  should  be 
regular  “bar”  shape  (see  Fig.  6,  left,  top),  but  as  the  tilt  angle 
is  increased  there  is  a  clear  bowing  effect  in  the  axial  direction 
which  corroborates  the  results  in  [7].  The  disks  are  almost 
unrecognizable  in  the  30-60°  tilts. 

Fig.  7  compares  the  results  of  the  mini-Defrise  phantom  in¬ 
serted  into  the  Tillable  breast  phantom.  Results  from  the  larger 
Defrise  phantom  indicated  extreme  axial  blurring  for  tilt  angles 
greater  than  15°.  This  experiment  investigated  smaller  tilt  angles 
to  evaluate  the  effects  of  incomplete  sampling  with  small  tilt  an¬ 
gles,  and  also  to  attempt  to  minimize  the  blurring  effect  but  still 
maximize  chest  wall  imaging  without  additional  orbits  to  com¬ 
plete  the  data.  A  bowing  effect  similar  to  that  seen  in  Fig.  6  is 
clear  in  these  reconstructions,  although  not  as  pronounced. 

Profiles  drawn  axially  across  the  left-center  of  the  breast  (see 
Fig.  8)  indicate  an  increase  in  counts  toward  the  nipple  due  to 
the  attenuating  effects  of  the  Defrise  phantom  within  the  breast. 


The  profiles  also  show  a  definite  axial  “stretching”  effect,  as 
the  15°  image  is  roughly  0.5  cm  longer  than  the  0°  image.  This 
result  indicates  that  visualization  of  axially  juxtaposed  tumors 
or  shell-like  tumors  may  be  difficult  to  discern. 

B.  Lesion  Imaging  with  Various  Tilt  Angles 

Reconstructed  sagittal  views  of  the  Tillable  breast  phantom 
containing  two  hot  lesions  combined  with  the  anthropomor¬ 
phic  torso  phantom  are  presented  in  Fig.  9  with  corresponding 
profiles  at  the  indicated  locations.  Background  heart,  liver,  and 
torso  contamination  were  present,  and  all  profiles  were  drawn 
across  the  center  of  the  lesions.  Table  I  compares  the  SNR  and 
contrast  values  for  each  lesion  and  tilt  angle.  Fig.  9  clearly  shows 
an  increase  in  both  contrast  and  SNR  as  tilt  angle  is  increased 
from  0°  to  30°.  Tilt  angles  beyond  30°  do  not  significantly  im¬ 
prove  either  contrast  or  SNR  values. 

The  third  lesion,  attached  to  the  anthropomorphic  torso 
phantom  anterior  to  the  chest  wall  and  superior  to  the  left 
breast,  was  also  reconstructed.  Although  not  shown,  the  lesion 
is  clearly  visible  in  all  reconstructions. 
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Fig.  7.  Single  slice  sagittal  view  OS-EM  reconstructions  of  mini-Defrise 
phantom  combined  with  Tillable  breast  phantom  at  indicated  angles.  Profile 
markers  are  drawn  axially  across  left-center  of  breast  phantom. 


Fig.  8.  Profiles  drawn  across  single  slice  sagittal  views  of  mini-Defrise 
phantom  (see  Fig.  7).  Left  corresponds  to  bottom  portion  of  the  images  shown 
in  Fig.  7.  Camera  head  was  tilted  away  from  chest  wall.  Note  increase  in 
breast’s  axial  length  with  greater  tilt  angle. 


IV.  Discussion  and  Conclusion 

Because  of  insufficient  sampling,  axial  blurring  is  evident  in 
all  TH-SPECT  reconstructions.  As  shown  with  the  larger  De- 
frise  phantom,  large  tilt  angles  (45°  and  60°)  yield  extremely 
distorted  images.  Because  the  ROR  can  be  decreased  as  tilt 
angle  is  increased,  it  is  important  to  find  an  acceptable  combi¬ 
nation  between  the  two,  specifically  one  which  minimizes  axial 
blurring  effects  while  still  maximizing  SNR  and  contrast  values. 
Tradeoffs  between  higher  resolution  parallel  hole  collimators 
versus  higher  sensitivity  collimators  should  be  considered.  The 
use  of  nonparallel  hole  collimators  should  also  be  considered 
in  this  application  since  these  geometries  could  potentially  help 
complete  undersampled  regions  of  the  Orlov  sphere  [8], 


Center  Lesion  Offset  Lesion 


Planar  Planar 


Fig.  9.  Sagittal  view  reconstructions  of  Tillable  breast  phantom  attached 
to  anthropomorphic  torso  phantom  containing  two  hot  lesions  (one  lesion 
axially  centered,  one  lesion  anterior  to  the  chest  wall  and  positioned  toward 
the  medial-superior  side  of  the  breast)  at  0°,  15°,  30° ,  and  45°  tilt  angles  as 
well  as  conventional  high-count  planar  images.  Profiles  for  each  lesion  are  also 
shown.  The  cold  spot  directly  behind  the  central  lesion  is  the  acrylic  support 
used  to  hold  lesion  in  place  and  the  cold  gap  between  the  breast  and  torso 
represents  phantom  chest  wall. 


As  Fig.  9  shows,  larger  tilt  angles  can  “stretch”  the  breast 
axially.  This  stretching  phenomenon  also  affects  the  lesion, 
and  can  potentially  blur  together  axially  offset  lesions  or 
distort  shell-like  lesions.  This  phenomenon  has  very  little 
effect,  however,  on  the  15°  degree  reconstruction.  Compared 
to  the  0°  reconstruction  which  has  no  distortions,  SNR  and 
contrast  values  between  the  two  are  significantly  different, 
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TABLE  I 

SNR'  and  Contrast  Values  for  Center  Lesion  and  Offset  Lesion  at 
0°,  15°,  30°,  and  45°  Tilt  Angles  as  well  as  a  Conventional 
High-Count  Planar  Image  (Craniocaudal  View) 


SNR 

[(Plesini!  “  pfrfwtt  )/  CJtmsisj] 

Contrast 

[( Plosion  “  |ibrcas»y  pbfwtst] 

ROR 

(cm) 

Center  lesion 

Planar 

2.10 

0.20 

18.1 

0° 

2.09 

0.25 

18.1 

15° 

3.56 

0.49 

10.5 

30° 

4.89 

0.59 

9.4 

45° 

5.15 

0.59 

5.9 

Offset  lesion 

Planar 

2.04 

0.20 

18.1 

0° 

2.13 

0.24 

18.1 

15° 

3.74 

0.51 

10.5 

30° 

6.32 

0.63 

9.4 

45° 

5.16 

0.49 

5.9 

providing  lesion  SNR  and  contrast  improvements  of  up  to  two 
times  for  the  center  lesion  (Table  I).  Although  30°  and  45° 
tilt  reconstructions  have  even  better  SNR  and  contrast  values, 
distortions  of  the  breast  are  much  more  evident,  compared  to 
the  smaller  tilt  angles. 

As  Table  I  indicates,  effective  radius  of  rotation  significantly 
decreases  as  tilt  angle  is  increased.  Although  the  45°  radius 
of  rotation  was  almost  half  that  of  the  30°,  extreme  insuffi¬ 
cient  sampling  as  well  as  attenuation  prevented  further  contrast 
and  SNR  improvements.  Less  dramatic  decreases  in  ROR  due 
to  increasing  tilt  angle  would  be  achieved  using  a  small-area 
dedicated  breast  SPECT  system  [5],  although  results  indicate 
SNR  and  contrast  values  should  improve  with  any  reduction  of 
ROR.  In  addition  to  viewing  lesions  near  the  torso-breast  in¬ 
terface,  TH-SPECT  reconstructions  may  considerably  improve 
SNR  and  contrast  values  even  for  lesions  within  the  breast,  pre¬ 
sumably  by  allowing  for  a  reduced  ROR.  Localization  of  the 
lesion,  however,  may  be  difficult  with  the  present  artifacts. 

The  blurring  and  distortion  artifacts  become  more  apparent 
with  greater  axial  nonuniformity  in  the  breast,  as  shown  in 
Fig.  7  with  the  mini-Defrise  phantom.  The  15°  tilt  reconstruc¬ 
tion  of  the  mini-Defrise  phantom  inside  the  tillable  breast 
phantom  highlights  a  definite  bowing  effect.  These  results 
indicate  that  artifacts  will  be  more  present  when  objects  are 
axially  aligned  and  imaged  with  a  TH-SPECT  technique.  This 
bowing  elfect  is  not  visualized  in  Fig.  9  reconstructions  simply 
due  to  the  uniformity  of  breast  tissue  and  spherical  lesions. 
Nonuniform  breast  lesions,  such  as  distorted  or  shell-like 
tumors,  may  result  in  further  image  degradation  even  at  small 
angular  tilts.  Although  dedicated  breast  SPECT  systems  under 
development  may  also  suffer  from  the  artifacts  illustrated  in 
this  study,  those  systems  may  further  provide  SNR  and  contrast 
improvements  due  to  reduced  RORs  [5]. 

When  the  reconstruction  was  limited  to  the  breast  only,  es¬ 
sentially  eliminating  background  interference  from  the  FOV,  a 
cupping  artifact  [9]  was  present.  This  cupping  artifact  is  an  ar- 


AOR 


Fig.  10,  Transverse  view  of  phantom  illustrating  TH-SPECT  acquisition  as 
well  as  an  additional  HAOR  orbit  to  gather  insufficiently  sampled  data. 

tificial  increase  in  activity  near  the  edge  of  the  FOV.  By  recon¬ 
structing  the  entire  image  space  (heart  and  liver  interference  in¬ 
cluded),  this  artifact  was  eliminated  [9].  It  remains  to  be  seen 
if  small-area  dedicated  breast  SPECT  systems  will  suffer  from 
these  limitations  and  if  reconstructed  images  from  a  small-area 
camera  contain  the  cupping  artifact. 

Possible  modifications  can  be  made  to  the  TH-SPECT  ge¬ 
ometry  to  account  for  insufficient  sampling  and  to  reduce,  if 
not  eliminate,  axial  blurring  artifacts.  Because  Orlov  Sphere  re¬ 
quirements  are  not  met  with  any  degree  of  camera  tilt  [3],  an 
additional  orbit  is  required  to  satisfy  these  conditions  for  fully 
quantitative  SPECT  imaging.  Fig.  10  introduces  one  such  ad¬ 
ditional  orbit.  If  TH-SPECT  VAOR  data  is  collected  and  com¬ 
bined  with  an  incomplete  HAOR  orbit  [10],  Orlov  Sphere  re¬ 
quirements  should  be  met  and  blurring  due  to  insufficient  axial 
sampling  should  be  eliminated.  This  short  arc-orbit  (ranging 
from  10°  to  90°  depending  upon  tilt  angle)  needs  only  to  gather 
enough  data  to  fully  sample  the  breast,  as  accurate  heart  and 
liver  reconstruction  may  not  be  necessary. 
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Abstract — Parallel-beam  tilted-head  single  photon  emission 
computed  tomography  (TH-SPECT)  was  previously  implemented 
on  a  SPECT  system  for  its  potential  to  image  breast  lesions 
and  nearby  axilla  of  seated,  upright  women.  All  TH-SPECT 
reconstructions  will  contain  artifacts  since  the  tilted  orbit  does  not 
satisfy  the  Orlov  sampling  criteria.  However,  it  is  not  clear  which 
reconstruction  method,  if  any,  is  better  suited  for  TH-SPECT 
data.  Here  a  geometric  derivation  of  the  ramp  filter  for  tilted 
parallel-beam  geometries  is  presented.  A  filtered  backprojection 
(FBP)  algorithm,  using  this  filter,  was  then  implemented  and  com¬ 
pared  with  an  iterative  ordered  subsets  expectation  maximization 
(OSEM)  algorithm,  using  TH-SPECT  data.  A  breast  scan  at  var¬ 
ious  tilt  angles  was  simulated  and  used  to  generate  a  noise  versus 
bias  study  for  both  methods.  Contrast  and  signal-to-noise  ratio 
(SNR)  values,  as  well  as  axial  elongation  present  in  all  TH-SPECT 
reconstructions,  were  characterized  using  a  mini-Defrise  disk 
phantom  placed  inside  a  tillable  breast  phantom  and  imaged  from 
0  to  15°  head  tilt.  A  tillable  breast  phantom  containing  lesions 
was  also  imaged  with  a  system  dedicated  to  prone  breast  SPECT 
from  0  to  30°  to  evaluate  the  effects  of  incomplete  sampling  due 
to  greater  tilt  angles.  FBP  noise  versus  bias  studies  indicated 
a  greater  increase  in  bias  with  tilt  angle  compared  to  OSEM 
reconstructions.  At  small  tilt  angles  about  the  mini-Defrise  disk 
phantom,  poorer  contrasts  were  obtained  with  FBP  compared  to 
OSEM  at  similar  noise  levels.  All  reconstructions  of  the  Tillable 
breast  phantom  indicated  axial  elongation  at  greater  tilt  angles, 
although  FBP  reconstructions  displayed  an  increase  in  stretching 
distortions  of  the  breast.  OSEM  SNR  and  contrast  values  were 
higher  at  all  degrees  of  tilt.  In  conclusion,  measured  results  indi¬ 
cate  OSEM  TH-SPECT  reconstruction  provides  better  contrast 
and  SNR  values  and  may  offer  better  shape  and  uniform  activity 
distribution  of  the  breast  compared  to  FBP  methods. 

Index  Terms — Breast  imaging,  filtered  backprojection  (FBP), 
ordered  subsets  expectation  maximization  (OSEM),  Orlov,  single 
photon  emission  computed  tomography  (SPECT). 


I.  Introduction 

TILTED-HEAD  single  photon  emission  computed  to¬ 
mography  (TH-SPECT)  parallel-beam  imaging  of  breast 
lesions  has  been  shown  to  offer  superior  lesion  contrast  and 
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signal-to-noise  ratio  (SNR)  values  compared  to  conventional 
SPECT  geometries  [1].  Reconstruction  of  TH-SPECT  data 
requires  a  three-dimensional  reconstruction  algorithm.  Iterative 
methods  have  thus  far  been  used  to  accurately  model  the 
complex  (albeit  incomplete)  sampling  of  frequency  space 
associated  with  the  tilted  circular  orbit.  It  may  be  useful, 
however,  to  reconstruct  TH-SPECT  data  using  analytic  filtered 
backprojection  (FBP)  methods  given  the  speed  and  linearity 
associated  with  such  methods. 

Conventional  two-dimensional  (2-D)  FBP  SPECT  recon¬ 
struction  for  parallel-beam  collimation  relies  on  the  standard 
one-dimensional  ramp  filter.  This  ramp  filter  has  also  been 
shown  to  be  the  backprojection  filter  for  three-dimensional 
(3-D)  tilted  parallel-beam  geometries  [2].  Due  to  insufficient 
Orlov  Sphere  sampling  [3],  TH-SPECT  acquisitions  introduce 
a  null  space  of  unmeasured  frequency  components  whose 
frequency-space  volume  is  dependent  upon  tilt  angle.  These 
unmeasured  frequency  components  may  have  negative  effects 
on  image  reconstruction.  One  such  effect  is  axial  stretching 
which  produces  unwanted  distortion  of  the  patient’s  body  con¬ 
tour  [1].  Another  related  effect  is  inaccurate  activity  distribution 
specifically  due  to  insufficient  sampling  of  the  Orlov  Sphere. 
With  these  negative  effects  in  mind,  it  is  unclear  which  image 
reconstruction  method  is  better  suited  for  TH-SPECT  data. 
This  paper  develops  a  geometric  derivation  of  the  ramp  filter 
for  TH-SPECT  and  then  compares  parallel-beam  TH-SPECT 
FBP  reconstruction  with  an  existing  iterative  ordered  subsets 
expectation  maximization  (OSEM)  algorithm  [1],  [4],  Simu¬ 
lation  studies  as  well  as  phantom  studies  are  performed  and 
evaluated. 


II.  Methods 

A.  TH-SPECT  Reconstruction 

Three  spatial  domain  coordinate  systems  are  used  (Fig.  1)  and 
are  defined  as  (.r,  y,  z),  (x\  yf,  z'\  and  (xf/,  y ",  zn).  The  ob¬ 
ject  is  expressed  within  the  fixed  (.r,  ?/,  z)  coordinate  system 
as  f(x,  y,  z).  The  (x/,  y\  zf)  system  is  a  rotational  coordinate 
system  rotated  about  the  z  axis  by  angle  0.  The  third  coordinate 
system  (.r",  yf\  z,r)  represents  the  tilt  of  the  ( x! ,  yr ,  zf)  system 
about  the  axis  by  tilt  angle  <j>.  The  object’s  frequency  domain 
representation  is  expressed  within  the  fixed  (//,  zz,  r)  system  as 
v,  r),  and  rotations  within  this  system  are  similarly  de¬ 
fined  as  (//,  i/,  r'),  and  (/z",  z/',  r"). 
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Fig.  1.  Illustration  of  tilted  geometry  and  (;r,  y.  c)  coordinate  system. 
Lower-right  illustration  indicates  breast  positioning  within  the  tilted  geometry, 
where  the  z  axis  passes  through  the  nipple  and  the  center  of  the  base  of  the 
breast.  Dashed  lines  indicate  axis  is  behind  object. 


Rotation  of  (.r,  y ,  z),  or  similarly  (/x,  v,  r),  by  angle  9  pro¬ 
duces  the  following  transfonnations: 

x!  =  x  cos  9  +  y  sin  0 

y  =  —x  sin  9  +  y  cos  0 

z'=z.  (1) 

Likewise,  rotation  of  (xf ,  y',  z'),  or  similarly  (pf,  v\  r'),  by 
angle  0  produces  the  following  transfonnations: 

xn  —  x  cos  9  cos  (j>-\-  y  sin  9  cos  cf>  —  z  sin  0 
yn  =  —x  sin  9  -f  y  cos  9 

z"  —  x  cos  9 sin  0  +  y  sin  9 sin  (j>  +  z  cos  0.  (2) 

1 )  Backprojection  Filter  for  N on-lilted  Parallel-Beam  Ge¬ 
ometry:  In  order  to  suggest  a  filter  for  TH-SPECT  data,  we  first 
examine  filtering  for  nontilted  parallel-beam  geometries.  The 
Fourier  projection-slice  theorem  implies  that  by  taking  projec¬ 
tions  at  multiple  angles  9  with  0  =  0,  values  of  F(/z,  v,  r)  will 
be  known  along  radial  planes  centered  at  the  origin  and  perpen¬ 
dicular  to  the  projections  [5]-[10].  Fig.  2  (top)  illustrates  the 
frequency  domain  representation  of  a  single  nontilted  projec¬ 
tion  at  angle  9 ,  with  u)  being  the  distance  from  one  arbitrary 
sample  point  to  the  axis  of  rotation  (AOR)  r.  Rotation  about  r 
fills  F(/x,  7v,  r)  with  multiple  circles  of  discrete  sample  points. 
Fig.  2  (bottom)  illustrates  a  slice  at  a  constant  value  of  r.  As 
the  radial  distance  lj  increases,  the  density  of  sample  points 
within  a  given  rotation  circle  decreases.  Specifically,  the  density 
of  sample  points  within  each  circle  is  directly  proportional  to 
1  fu9  the  inverse  of  the  circles’  circumference,  or,  in  our  rotated 
coordinate  system,  1  /z/7.  Because  the  distance  7  between  all  ro¬ 
tation  circles  is  constant,  the  density  of  circles  within  the  plane 
is  also  constant.  The  density  of  sample  points,  which  can  be  ex¬ 
pressed  as  the  density  of  circles  times  the  density  of  points  per 
circle,  is  therefore  proportional  to  1  /lj.  The  filter  to,  which  is  the 
well-known  ramp  filter,  compensates  for  this  varying  sampling 
density.  Because  lower  frequencies  are  sampled  more  densely 
than  higher  frequencies,  a  filter  is  needed  to  correct  for  this  dif¬ 
ference  in  sampling.  Because  each  slice  along  r  is  identical 
in  sample  density,  the  backprojection  filter  for  nontilted  par¬ 
allel-beam  data  is  one-dimensional  across  projections. 


Axis  of  rotation 


Fig.  2.  (Top)  Illustration  of  Fourier  projection-slice  theorem  for  nontilted 
data.  The  plane  indicates  the  frequencies  sampled  by  a  projection  at  angle  9. 
cj  represents  the  distance  between  an  arbitrary  sample  point  and  the  r  axis. 
(Bottom)  Frequency  domain  slice  at  a  constant  value  of  r  showing  the  discrete 
sample  points  at  distances  v"  —  uf  =  a;  from  the  AOR  and  rotated  by  9. 
Distance  7  represents  the  constant  distance  between  circles  of  discrete  sample 
points. 


2)  Backprojection  Filter  for  Tilted  Parallel-Beam  Geom- 
etiy:  For  tilted  parallel-beam  data  (0  ^  0),  Fig.  3  (top) 
illustrates  the  region  of  frequency  space  sampled  by  a  projec¬ 
tion  at  angles  9  and  0.  The  distance  w  from  the  axis  of  rotation 

o  '  ft - - - - Ty — — 

r  to  a  sample  point  is  u>  =  y  vn~  4-  r"  sin“  0. 

Sampled  frequencies  correspond  to  discrete  values  of  i/',  and 
each  value  of  vn  lies  on  a  circle  of  sample  points  having  a  radius 
u;,  with  one  sample  point  per  angle  9.  The  density  of  points  along 
a  given  circle  is  inversely  proportional  to  the  circles’  radius,  that 
is,  proportional  to 


yfi;"2  +  r"2  sin2  0 

In  contrast  to  nontilted  data,  the  distance  7  between  circles 
in  this  plane  is  not  constant  but  rather  depends  upon  v"  and  rn 
(Fig.  3,  bottom).  The  density  of  circles  is  proportional  to  the 
reciprocal  of  the  distance  between  circles,  that  is,  proportional 
to 

dvn  yj uff2  -f  t'/2  sin2  0 

__  =  __ 


(4) 
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Fig.  3.  (Top)  Illustration  of  Fourier  projection-slice  theorem  for  tilted  data. 
Plane  illustrates  the  frequency  space  sampling  for  a  projection  at  angles  0  and  o. 
u>  represents  distance  between  arbitrary  sample  point  and  the  r  axis.  (Bottom) 
Frequency  domain  slice  at  constant  value  r  ±  0  built  up  from  discrete  sample 
points  along  v”  and  over  rotation  angles  9.  Frequencies  are  not  measured  within 
the  concentric,  inner  circle  of  radius  of  t"  sin  if).  Unmeasured  frequency  region 
decreases  as  r  — +  0  and  becomes  Fig.  2  (bottom)  when  r  =  0.  Note  distance 
is  now  a  function  of  both  v”  and  r"  sin  (f>.  Distance  7  represents  the  distance 
between  circles  built  up  from  discrete  sample  points.  Note  circles  are  more  dense 
near  u)  =  r"  sin  (f>  and  less  dense  with  increasing  u,*.  For  simplicity,  sample 
points  were  only  shown  at  45-degree  increments.  Conventional,  finer  sampling 
would  diminish  this  nonuniform  density  appearance  within  each  circle. 


As  vn  increases,  the  distance  between  neighboring  circles  in¬ 
creases.  These  circles  are  more  dense  near  u  =  r,f  sin  (j>  than 
at  higher  gj.  For  uj  <  | r"  sin  <j>\,  the  density  of  sample  points 
within  the  plane  is  zero.  For  w  >  |r"  sin  <j> |,  the  density  of 
sample  points  within  the  plane  is  the  product  of  (3)  and  (4), 
which  is  1 

Finally,  because  the  spacing  of  sample  points  in  the  r  direc¬ 
tion  for  tilted  data  is  constant,  although  compressed,  a  third  con¬ 
stant  component,  1  /  cos  </>,  can  be  included. 

Multiplying  these  three  density  components  yields 

1  yj i/"2  +  t"2  sin"  1  _  1 

Jv"'1  +  t"2  sin2  <b  Kl  cos  <{>  | i/rr  \  cos  <j> 

(5) 

for  the  density  of  all  points  when  u  >  \r"  sin  <f>\.  This  density 
is  the  modulation  transfer  function  (MTF)  of  the  ideal  (no  atten¬ 
uation,  perfect  spatial  resolution)  tilted  parallel-beam  measure¬ 
ment.  Fora;  <  |r"  sin  <f>\,  the  MTF  is  zero,  so  those  frequencies 


Fig.  4.  Illustration  of  missing  conical  regions  present  in  all  TH-SPECT 
frequency  space  representations.  Volume  of  missing  region  increases  with  tilt 
angle  0. 


Area  imaged 

Fig.  5.  Transverse  view  of  digital  anthropomorphic  Peter  phantom.  Squared 
breast  region  indicates  the  subsection  of  phantom  which  was  imaged. 

cannot  be  recovered.  These  frequencies  can  be  represented  as 
missing  conical  regions  whose  volumes  increase  with  tilt  angle 
0  (Fig.  4). 

For  a;  >  \rn  sin  <f> |,  (5)  indicates  that  the  MTF  can  be  in¬ 
verted  by  the  scaled  ramp  filter,  \vn  cos  <j> |.  This  geometrical 
derivation  is  consistent  with  previous  derivations  of  the  back- 
projection  filter  for  the  tilted  parallel-beam  geometry  [2]. 

B.  Simulated  Breast  Scan 

In  order  to  compare  OSEM  and  FBP  for  tilted  parallel-beam 
breast  imaging,  a  vertical  axis  of  rotation  (VAOR)  scan  was 
simulated  using  a  section  of  the  digital  anthropomorphic  Peter 
phantom  [11].  This  computer  phantom  consists  of  various 
volumetric  regions  in  the  torso  (Fig.  5).  For  this  study,  only 
the  left  breast  (11  cm  base-to-nipple,  14  cm  base  diameter) 
containing  a  uniform  activity  concentration  was  scanned  and 
no  torso  activity  was  present.  The  phantom  was  implemented 
on  a  168  x  168  x  168  grid,  voxels  0.35  cm  on  a  side.  The 
camera’s  AOR  coincided  with  a  line  from  the  base  of  the  breast 
through  the  nipple  (Fig.  6).  The  simulated  VAOR  scan  used 
a  22.4  x  22.4  cm  field-of-view  (FOV)  parallel-hole  gamma 
camera,  and  each  projection  data  set  was  acquired  for  128 
angular  views  over  360°.  Each  projection  was  64  x  64  pixels 
and  pixel  size  was  0.35  cm2.  Projection  data  were  simulated 
for  0°,  15°,  30°,  and  45°  tilt  angles.  Expected  values  for 
VAOR  projection  data  were  analytically  computed  based  on 
forward  projections  of  activity  within  the  breast  region  only. 
Compton  scattering,  attenuation,  and  collimator  and  crystal 
spatial  resolution  were  not  modeled. 
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Camera  head 


Fig.  6.  Illustration  of  simulated  breast  scan.  Breast  represents  subsection  of 
digital  anthropomorphic  Peter  phantom  and  O  represents  angular  camera  tilt. 


Fig.  7.  Schematic  of  breast  phantom  with  three  mini-Defrise  disks.  Nipple  is 
at  top  and  curvature  at  bottom  corresponds  to  torso  shape. 


Fig.  8.  Schematic  of  950  ml  breast  phantom  containing  lesions.  Large  lesion 
is  located  at  a  superior  medial  aspect  of  the  breast,  6  cm  from  nipple  and  3  cm 
from  a  normal  to  the  surface  of  the  breast.  Smaller  lesion  is  located  at  an  inferior 
medial  aspect  of  the  breast,  6  cm  from  nipple  and  3  cm  from  a  normal  surface  of 
the  breast  (top  schematic  shows  relative  height  and  lesion  displacement;  Bottom 
schematic  illustrates  actual  lesion  orientation  in  the  breast). 

An  ensemble  study  was  generated  consisting  of  50  acqui¬ 
sitions  per  tilt  angle  and  based  on  Poisson  variation  about 
the  expected  values.  OSEM  reconstructions  were  performed 
with  the  tilted  geometry  modeled  in  the  photon  detection 
probabilities,  and  noise  and  bias  values  were  computed  for 
1-50  iterations  (eight  subsets).  Reconstructed  OSEM  images 
were  32  x  32  x  32  cubic  voxels  each,  0.7-cm  wide.  Compton 
scattering,  attenuation,  and  collimator  and  crystal  spatial 
resolution  were  not  modeled  within  the  reconstruction.  FBP 
reconstruction  was  similarly  performed  with  images  also 
32  x  32  x  32  cubic  voxels,  0.7-cm  wide.  FBP  images  were 
post-filtered  using  a  3-D  Butterworth  filter.  Butterworth  filter 
cut-off  frequencies,  the  frequencies  at  which  the  gain  equals 
0.707,  ranged  from  0.1-1. 0  cycles/pixel  with  frequency  order 
equal  to  10. 


Fig.  9.  Photograph  of  Tillable  breast  phantom  suspended  over  the  ASET 
system. 


Noise 

Fig.  10.  Noise  versus  bias  plots  for  OSEM  (top)  and  FBP  (bottom) 
reconstructions  of  simulated  breast  scan  at  0°,  15°,  30°,  and  45°  tilt.  OSEM 
data  points  correspond  to  iterations  1-50  in  increments  of  one  iteration  from 
left  to  right.  FBP  data  points  indicate  Butterworth  filter  cutoffs  (1.0  to  0.1 
cycles/pixel  in  increments  of  0.1  cycles/pixel  from  right  to  left,  where  furthest 
right  point  indicates  no  post-filtering).  OSEM  yielded  the  lower  bias  for  similar 
noise  values. 

Noise  versus  bias  curves  were  generated  for  both  OSEM  and 
FBP  reconstructions.  Noise  was  computed  as  the  standard  devi¬ 
ation  over  the  ensemble  on  a  pixel-by-pixel  basis.  Bias  was  de¬ 
fined  as  J  Eili  ~  arf*  )2 ,  where  N  is  the  number  of  voxels, 

Xi  is  the  it h  pixel  value  of  the  average  image,  and  x\  is  the 
value  of  the  ith  pixel  in  the  true  phantom  image. 

C.  Defrise  Disk  Phantom  Scan 

Mini-Defrise  phantom  disks  (Data  Spectrum  Corp.,  Hills¬ 
borough,  NC)  were  placed  inside  a  Tillable  breast  phantom 
(Model  ECT/FL  BR/A,  Data  Spectrum  Corp.,  Hillsborough, 
NC)  (Fig.  7).  Data  were  acquired  using  one  head  of  a  tiltable 
two  headed  SPECT  system  (Vision  T-22,  Summit  Nuclear, 
Twinsburg,  OH)  and  a  parallel  beam  low  energy  high  resolution 
(LEHR)  collimator  (Hitachi,  Tokyo,  Japan)  (hole  length  — 
41.0  mm,  flat  to  flat  hexagonal  hole  size  =  1.8  mm,  septal 
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OSEM 


FBP 


Phantom 


Transaxial  Coronal 

15"  30°  45*  0°  15*  30^  45”  0*  15*  30°  45° 


Fig.  1 1 .  Transaxial,  coronal,  and  sagittal  (with  respect  to  the  breast  reference  frame)  OSEM  and  FBP  reconstruction  slices  for  simulated  breast  scan.  Slices  were 
taken  for  0°,  15°,  30°,  and  45°  tilt  and  corresponding  truth  phantom  slices  located  at  bottom. 


thickness  =  0.18  mm).  The  breast  phantom  contained  three 
cold  Defrise-phantom  disks  measuring  1 .0  cm  in  thickness  and 
8.0  cm  in  diameter  and  with  a  1.0-cm  gap  between  disks.  The 
breast  phantom  was  filled  with  water  having  a  22.5  /xCi/ml 
concentration  of  Tc-99m.  The  axis  of  the  breast  phantom  was 
then  aligned  parallel  to  the  camera’s  AOR  through  the  nipple, 
and  the  radius  of  rotation  (ROR)  was  measured  perpendicular 
to  the  center  of  the  camera  and  to  the  nipple.  Projection  data 
were  acquired  for  tilt  angles  of  0°,  5°,  10°,  and  15°  having 
RORs  of  12.1  cm,  1 1.8  cm,  1 1.4  cm,  and  10.5  cm,  respectively. 
The  ROR  decreased  with  greater  tilt  angle  due  to  greater  access 
to  the  breast.  Total  acquisition  time  was  10  s  per  angular  step 
or  21.33  min  over  360°.  Acquisition  times  were  adjusted  to 
compensate  for  radioactive  decay. 

OSEM  (eight  subsets,  three  iterations)  and  FBP  (Butterworth 
filtered,  cut-off  =  0.5  cycles/cm,  order  =  10)  reconstructions 
were  performed  for  all  tilt  angles,  and  final  images  were 
128  x  128  x  128  cubic  voxels,  0.48  cm  on  a  side.  Iteration  and 
filter  parameters  were  determined  from  aforementioned  noise 
versus  bias  study. 

D.  Fillable  Breast  Phantom  With  Lesions 

The  fillable  breast  phantom  (950  ml)  containing  two  spher¬ 
ical  lesions,  1.0  cm  diameter  (0.52  ml)  and  0.6  cm  diameter 
(0.11  ml)  (models  ECT/HOL-468/A,  Data  Spectrum  Corp., 
Hillsborough,  NC),  was  also  imaged  in  this  experiment  (Fig.  8). 
Data  were  acquired  using  a  <13  x  13  cm  FOV  LumaGEM™ 
gamma  camera  (Gamma  Medica,  Inc.,  Northridge,  CA)  (hole 
length  =  23.6  mm,  flat  to  flat  hexagonal  hole  size  =  1.22  mm, 
septal  thickness  =  0.2  mm)  as  a  component  on  a  dedicated 
Application  Specific  Emission  Tomograph  (ASET)  gantry  [12] 
(Fig.  9).  The  breast  phantom  contained  "mTc-pertechnetate 
(140.6  keV)  activity  and  water  and  lesions :  breast  activity 
concentration  was  5.94:1.  Images  were  acquired  for  128 
angular  views  over  360°,  and  each  projection  was  54  x  54 
pixels,  each  2.2  mm2.  Total  acquisition  time  was  20  min,  or 
9.375  s  per  angular  step.  Projection  data  were  acquired  for  0°, 
10°,  20°,  and  30°  tilt  angles  with  RORs  of  7.47  cm,  6.63  cm, 
5.48  cm,  and  5.48  cm,  respectively,  measured  from  the  camera 
head  to  the  center-of-rotation  point.  Greater  tilted  angles  allow 
for  the  smaller  RORs. 

OSEM  (eight  subsets,  three  iterations)  and  FBP  (Butterworth 
filtered,  cut-off  =  0.5  cycles/cm,  order  =  10)  reconstructions 


were  performed  for  all  tilt  angles,  and  final  images  were 
54  x  54  x  54  cubic  voxels,  0.22  cm  on  a  side.  Iteration  and 
filter  parameters  were  determined  from  aforementioned  noise 
versus  bias  study  to  provide  similar  noise  levels  for  both 
methods. 

For  both  reconstruction  methods  at  all  degrees  of  tilt,  an 
elliptical,  1 1  pixel  region-of-interest  (ROI)  was  taken  for  the 
large  lesion  and  a  similar  121-pixel  ROI  was  taken  in  the  breast. 
Lesion  Contrast  was  Calculated  as  (^lesion  ~~  Pbreast)/ ^breast,  y 
and  SNR  as  {fusion  ~  Abreast) /(?br east,  where  IX  =  mean  and 
a  =  standard  deviation  of  selected  ROIs. 

III.  Results 

A.  Simulated  Breast  Scan 

Fig.  10  presents  noise  versus  bias  plots  for  OSEM  and  FBP 
reconstructions  of  the  simulated  breast  scan.  Fig.  1 1  presents 
transaxial,  sagittal,  and  coronal  mid-slices  (with  respect  to  the 
breast  reference  frame)  for  FBP  (no  post-filtering),  and  OSEM 
(3  iterations)  reconstruction  methods  at  0°,  15°,  30°,  and  45°. 
Iteration  cut-offs  were  chosen  to  provide  approximately  equal 
noise  for  both  reconstruction  methods.  FBP  results  indicate  a 
greater  increase  in  backprojection  artifacts  with  greater  angular 
tilt.  OSEM  reconstructions  were  better  able  to  maintain  a  clear, 
defined  border  around  the  breast,  but  both  methods  displayed 
some  degree  of  axial  stretching  of  the  breast  as  well  as  inac¬ 
curate  activity  distribution  with  greater  tilt.  Stretching  effects 
were  not  found  transaxially,  however,  due  to  complete  frequency 
sampling  of  the  r  =  0  plane  (Fig.  3).  This  is  consistent  with 
tilted  cone-beam  reconstruction  [13],  [14],  where  distortions 
were  most  pronounced  along  the  direction  parallel  to  the  ro¬ 
tation  axis  but  not  along  the  transverse  direction.  Ideally,  the 
breast  would  be  fully  angularly  sampled  to  satisfy  the  Orlov 
condition,  and  hence  reduce  currently  observed  elongation  and 
object  distortion  artifacts  [1],  [4]. 

FBP  noise  versus  bias  plots  indicate  greater  bias  with  tilt 
angle  compared  to  OSEM  reconstructions  (Fig.  10).  A  50%  in¬ 
crease  in  bias  was  noted  between  FBP  0°  and  15°  recontruc- 
tions,  while  only  a  15%  increase  in  bias  was  noted  between 
OSEM  0°  and  15°  reconstructions,  compared  at  similar  noise 
levels.  Butterworth  filtering  tended  to  offer  the  best  compro¬ 
mise  between  noise  and  bias  for  FBP  reconstructions,  specifi¬ 
cally  at  a  0.3  cycles/pixel  cutoff.  Three  iterations  displayed  the 
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Fig  12  Sagittal  view  OSEM  and  FBP  reconstructions  of  mini-Defrise  Fig-  13.  Sagittal  view  OSEM  and  FBP  reconstructions  of  anthropomorphic 
phantom  combined  with  flllable  breast  phantom  and  corresponding  profiles.  torso  phantom  combined  with  tillable  breast  phantom  containing  lesions.  Profile 
Profile  markers  are  drawn  axially  across  left-center  of  breast  phantom.  markers  are  drawn  axially  across  left-center  of  breast  phantom. 


TABLE  I 

Contrast  and  SNR  Values  of  Large  Lesion 
for  OSEM  and  FBP  Reconstructions 


SNR 

Contrast 

ROR  (cm) 

OSEM 

0° 

9.98 

2.65 

7.47 

10° 

13.43 

3.84 

6.43 

20° 

12.23 

2.82 

5.48 

30° 

12.85 

3.57 

5.48 

FBP 


0° 

5.57 

2.26 

7.47 

10° 

6.39 

2.13 

6.43 

20° 

6.59 

2.75 

5.48 

30° 

9.10 

3.25 

5.48 

best  noise/bias  trade-off  for  OSEM  0°,  15°,  and  30°  data  and 
four  iterations  were  found  optimal  for  45°  data.  OSEM  recon¬ 
structions  displayed  the  lower  overall  bias. 


B.  Defrise  Disk  Phantom  Scan 

Fig.  12  displays  reconstructions  of  mini-Defrise  phantom 
disks  placed  inside  a  flllable  breast  phantom.  The  resulting 
profiles  indicate  a  bowing  effect  slightly  visible  at  10°  tilt,  and 
more  clearly  visible  at  15°  tilt.  OSEM  Defrise  disk  contrast 
was  significantly  better  at  greater  tilt  angles  compared  to 
FBP  results.  OSEM  reconstruction  sets  displayed  a  one  pixel 
(0.48  cm,  3.8%)  increase  in  breast  length  at  15°  tilt,  measured 
from  the  nipple  to  the  anterior  wall.  FBP  profiles  indicated  a 
two  pixel  (0.98  cm,  7.6%)  increase  in  breast  length.  All  profiles 
drawn  axially  across  the  right-center  of  the  breast  indicate  an 
increase  in  counts  toward  the  nipple  due  to  attenuating  effects 
of  the  Defrise  disk  phantom  within  the  breast. 

C.  Fillable  Breast  Phantom  With  Lesions 

Table  I  presents  contrast  and  SNR  values  for  reconstructions 
of  the  flllable  breast  phantom  containing  lesions  (Fig.  13). 
OSEM  reconstructions  showed  greater  SNR  and  contrast  values 
at  all  tilt  angles.  Because  camera  tilt  allows  greater  access  to  the 
breast  and  in  turn  decreases  the  effective  ROR,  overall  contrast 
and  SNR  values  generally  increased  with  tilt  angle,  confirming 
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the  effectiveness  which  can  be  achieved  with  TH-SPECT  and 
dedicated  compact  gamma  cameras. 

IV.  Discussion  and  Conclusion 

This  paper  has  presented  a  geometric  derivation  of  the  ramp 
filter  fortified  parallel-beam  geometries.  FBP  simulation  results 
indicated  a  greater  increase  in  artifacts  with  increasing  tilt,  as 
compared  with  iterative  reconstruction,  which  better  maintained 
overall  breast  shape.  Simulation-based  noise  versus  bias  studies 
indicated  near-optimal  iteration  and  cut-off  parameters  for  each 
method  and  showed  less  bias  at  a  given  noise  for  OSEM  as  com¬ 
pared  to  FBP.  Phantom  studies  demonstrated  consistently  higher 
SNR  and  contrast  values  for  OSEM  reconstructions  at  all  tilt 
angles. 

TH-SPECT  data  do  not  sufficiently  sample  the  Orlov  sphere 
and  therefore  artifacts  are  present  in  all  TH-SPECT  reconstruc¬ 
tions.  The  effect  of  these  insufficiencies  can  be  examined  in  fre¬ 
quency  space  as  missing  conical  regions  shown  in  Fig.  4.  This 
“missing  cone”  problem  has  been  investigated  in  a  variety  of 
different  3-D  imaging  systems  with  limited  view  angles  [2], 
[15]— [1 7].  If  frequency  components  within  the  missing  cone 
and  especially  along  the  r  axis  are  assumed  to  be  zero,  the  re¬ 
constructed  mean  densities  of  all  planes  perpendicular  to  £  will 
be  equal.  This  has  the  effect  of  producing  unwanted  activity  in 
areas  where  activity  should  be  zero  or  near-zero.  These  activity 
distortions  are  observed  in  FBP  techniques  as  streak  artifacts 
(Fig.  1 1)  along  the  axis  of  rotation.  Due  to  complete  sampling 
of  the  r  =  0  plane,  transverse  slices  do  not  exhibit  streak  arti¬ 
facts,  although  in-plane  distortions  are  present. 

The  quality  of  TH-SPECT  reconstructions  clearly  degrades 
with  tilt  angle,  specifically  due  to  the  volumetric  increase  in  null 
space.  The  fraction  of  null  space  to  total  frequency  volume  is 
given  as 

Null  space 

Total  volume,  ~  ^  (  ) 

With  30°  tilt,  only  13.4%  of  frequency  space  accounts  for 
missing  data.  With  60°  tilt,  however,  half  of  the  frequency  data  is 
missing.  This  nonlinear  increase  in  unsampled  volume  indicates 
tilt  angles  over  30°  may  not  be  useful  for  breast  SPECT. 

The  conjugate  views  for  nontilted  projections  exist,  and, 
therefore,  9  need  only  be  sampled  from  0  to  7 r.  For  tilted 
projections,  however,  the  conjugate  view  does  not  exist  and  9 
should  be  sampled  from  0  to  27r.  This  is  illustrated  in  Fig.  3 
(bottom)  as  rotation  from  0°  to  180°  degrees  does  not  even 
fully  sample  the  frequency  space  region  outside  the  circle  of 
radius  r"  sin  (j> . 

Herein,  we  have  investigated  tilted  orbits  in  which  the  camera 
center  moves  along  a  circle,  elucidating  the  nature  of  the  incom¬ 
plete  sampling  provided  by  these  orbits  and  the  effectiveness  of 


FBP  and  OSEM  in  utilizing  this  incomplete  data.  In  other  work, 
we  proposed  circle-plus-arc  [1],  [4]  and  more  elaborate  orbits 
that  do  provide  complete  sampling  of  an  isolated  breast,  and  we 
have  developed  a  SPECT  system  capable  of  implementing  such 
orbits  [12]. 
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The  objective  of  this  study  was  to  investigate  the  effect  on  contrast  and  SNR  of  shielding 
the  bed  with  lead.  Shielding  reduces  other-organ  photon  contamination,  but  may  result  in 
incomplete  sampling  of  other  organs  and  increased  data  inconsistency,  if  shielding 
attenuation  is  not  modeled  dining  reconstruction.  A  Trionix  XLT  scanner  with  one 
pinhole  collimator  was  used  to  obtain  experimental  data  on  a  breast  phantom  mounted  on 
a  torso  phantom.  The  torso  phantom  contained  the  liver,  lungs,  and  myocardium,  all  filled 
with  varying  amounts  of  99mTc.  The  breast  phantom  contained  two  spherical  lesions, 
also  filled  with  99mTc.  The  bed  contained  a  cut-out  for  the  prone  breast  attachment.  Ten 
scans  were  performed  with  bed  shielding,  and  without,  in  order  to  determine  the 
statistical  significance  of  the  analysis.  The  incomplete  circular  orbit  angular  range  was 
186°  and  consisted  of  94  views  with  a  2°  step  between  each  view.  The  data  were  then 
reconstructed  by  OSEM  with  8  subsets,  using  up  to  30  iterations  without  modeling  shield 
attenuation.  A  3-D  region-of-interest  analysis  was  performed  on  a  lesion  to  evaluate 
shielding  effects.  The  use  of  shielding  slightly  reduced  noise  in  contrast  estimates,  from 
about  0.25  to  0.23.  However,  shielding  increased  bias  in  contrast  from  0.24±0.08  to 
0.44±0.07.  SNRs  were  also  lower  with  shielding.  These  results  suggest  that  if  shielding  is 
to  be  of  benefit,  shielding  attenuation  must  be  modeled  within  iterative  reconstruction. 


Lead  Shielding  in  Pinhole  Breast  Imaging 

T.  Lee,  K.*Braun,  R.  Jaszczak,  J.  Bowsher,  K.  Bobkov 

Introduction: 

When  imaging  a  breast  with  a  large  field-of-view  camera, 
a  pinhole  collimator  is  advantageous  in  its  ability  to  decrease  the 
radius-of-rotation  (ROR)  and  get  close  to  the  breast.  It  has  been 
shown  that  SPECT  imaging  of  a  breast  phantom  using  a  pinhole 
collimator  has  significantly  reduced  the  effects  of  breast- tissue 
activity  that  overlies  and  underlies  the  lesion.  [1]  However,  when 
the  breast  is  mounted  on  an  anthropomorphic  torso  phantom, 
increased  photon  contamination  is  introduced  from  the 
myocardium,  liver  and  gall  bladder  that  may  impact  the  clarity  of 
the  breast  image.  The  use  of  a  pinhole  collimator  may  obtain  data 
with  decreased  photon  contamination  from  these  organs. 

Although  the  organ  of  interest  is  the  breast,  the 
myocardium,  liver,  and  gall  bladder  also  take  up  the  ""’Tc-MIBI. 
Thus,  in  our  phantom  studies,  we  are  investigating  the  difference 
that  lead  shielding  will  make  in  the  contribution  of  photons  from 
the  outlying  organs  that  are  detected  through  the  breast  phantom 
by  the  pinhole  collimator.  More  specifically,  we  are  interested  in 
determining  how  the  radioactivity  uptake  in  these  organs  affects 
the  contrast,  bias  and  SNR  in  the  phantom  breast  images  using 
region-of-interest  (ROI)  analysis. 

Data  Acquisition: 

A  Trionix  XLT  scanner  (Trionix  Research  Lab, 
Twinsburg,  Ohio)  with  one  pinhole  collimator  was  used  to  obtain 
the  experimental  data  on  a  breast  phantom  mounted  on  an 
anthropomorphic  torso  phantom  placed  in  a  prone  position  on  the 
bed.  Two  sets  of  scans  were  performed,  each  containing  10 
individual  scans.  The  first  set  of  scans  consisted  of  the  torso 
phantom  and  all  its  organs  filled  with  activity  as  shown  in  Figure  1. 
The  second  set  of  scans  consisted  of  the  same  configuration  of  the 
torso  phantom  as  the  first  scan,  except  with  the  addition  of  lead 
shielding  placed  over  the  entire  table  except  for  where  the  breast 
was  located. 


Figure  1 :  A)  ICO  angular  range  B)  Experimental  Setup 


These  scans  had  an  angular  range  from  106°  to  292°.  All 
the  scans  consisted  of  an  ROR  of  12.9  cm,  a  tungsten  pinhole 
aperture  of  3  mm,  and  pinhole  focal  length  of  20  cm.  All  the 
scans  consisted  of  94  views  with  2°  steps  between  each  view.  The 
scan  time  per  view  for  the  first  scan  was  3  seconds.  The  scan  time 
per  view  for  each  scan  thereafter  was  adjusted  in  length  for 
radioactive  decay  so  that  the  effective  scan  times  were  the  same 
for  all  the  scans.  A  total  of  20  scans  were  performed.  Ten  scans 
of  each  configuration  were  performed  to  obtain  better  statistical 
results.  The  projection  data  was  stored  in  a  grid  size  of  256x128 
(pixel  size  =  1 .72mm)  for  further  processing. 

Phantom  Study: 

The  torso  phantom  (Model  ECT/TOR/P,  Data  Spectrum 
Corp.,  Hillsborough,  NC)  contained  the  lungs,  liver,  and  heart. 

The  breast  was  mounted  on  the  left  side  of  the  torso  and  placed  in 
a  prone  position  on  the  imaging  bed.  Two  hollow  spheres,  used  to 
simulate  lesions,  were  placed  in  the  breast,  one  approximately 
three  centimeters  below  the  nipple  and  the  other  close  to  the  chest 


wall  towards  the  sternum.  Both  had  an  outside  diameter  of  2.2  cm. 
This  relatively  large  sphere  size  was  chosen  in  order  to  assist  with 
positioning  the  ROI’s  in  the  reconstructed  image.  The  liver 
(volume  =  1180  mL)  was  filled  with  22.5  pCi/mL,  the 
myocardium  (volume  =  118  mL)  was  filled  with  30  pCi/mL,  and 
the  lesions  (volume  =  5.5  mL)  were  filled  with  7.5  pCi/mL.  The 
remaining  parts  of  the  torso  and  breast,  the  ventricle,  and  the  lungs 
were  filled  with  1.5  pCi/mL.  To  expedite  the  scanning  time, 
increased  concentrations  above  normal  clinical  concentrations 
were  used. 


Data  Processing: 

In  this  initial  study,  the  projection  data  were 
reconstructed,  without  shielding  attenuation  modeling  or  scattering 
and  attenuation  correction,  with  an  OS-EM  (Ordered  Subset 
Expectation  Maximization)  algorithm  for  faster  processing  with 
minimal  image  quality  compensation.  Thirty  iterations  with  an  OS 
level  (or  the  number  of  ordered  subsets)  of  8  were  used  over  94 
projection  views.  The  reconstructed  image  was  in  a  grid  size  of 
128x128x100  with  a  voxel  size  of  2mm 

Results  and  Discussion: 

The  reconstructed  images  for  individual  scans  were 
relatively  noisy  making  the  ROI  delineation  difficult.  One  method 
to  overcome  this  difficulty  was  to  sum  up  all  the  individual  scans 
for  each  set  before  reconstruction  to  simulate  a  long  scan, 
assuming  that  the  phantom  did  not  move  during  data  acquisition. 
The  simulated  long  scan  provided  for  easier  ROI  delineation  since 
noise  was  reduced.  This  is  shown  in  Figure  2.  Comparison  of  the 
profiles  along  the  white  line  shows  that  the  individual  scan  is 
noisier  than  the  long  scan. 


a)  Individual  Scan  b)  Long  Scan 

Figure  2:1s*  Iteration  -  OSEM 


ROI  analysis: 

The  sphere  closer  to  the  chest  wall  was  distorted  in  shape 
after  reconstruction,  most  likely  due  to  insufficient  sampling.  The 
circled  lesion  as  shown  in  Figure  2  is  used  for  analysis.  Contrast 
and  SNR  (signal  to  noise  ratio)  values  were  employed  to  evaluate 
shielding  effects.  They  can  be  mathematically  expressed  as, 


ContrastJl*dtMl  • 


IAL-AT, 


where  NROi  and  NBG  are  the  mean  pixel  intensity  in  the  ROI  and 


the  background  respectively,  and  abgis  the  standard  deviation  in 


pixel  intensity  over  the  ROI. 

A  spherical  ROI  of  0.9  cm  in  radius  was  created  within 
the  lesion  volume  with  its  center  at  the  centroid  of  the  lesion. 


Statistical  evaluation  was  highly  dependent  upon  where  the 


background  was  because  scattering  and  attenuation  were  not 
corrected  in  the  image.  To  achieve  uniformity,  a  background 
mask  was  obtained  by  creating  a  shell  of  1 .2cm  in  thickness 
around  the  sphere  ROI  with  a  gap  (0.8  cm)  between  the  surface  of 
the  ROI  mask  and  the  inner  surface  of  the  background  shell.  The 


masks  were  applied  to  the  reconstructed  images  from  the  short 


> 

scans  for  statistical  evaluation  as  shown  in  Table  1  and  Figure  3. 
v‘Mean’  value  (for  both  contrast  and  SNR)  is  the  average  of  10 
different  scans.  ‘Noise’  value  is  an  estimate  of  the  standard  error 
of  the  distribution  from  which  contrast  or  SNR  rises.  Bias  is  the 
difference  between  the  measured  contrast  and  the  true  contrast 
(Qrue  =  4). 


Table  1 :  Shielding  Effect 


Without  Shielding  -  No  Filter  Applied  | 

Contrast 

_ l 

Iteration 

Noise 

Bias 

1th 

USES 

0.19 

0.7*0.06 

4.5*0.28 

E£B 

2nd 

3.8*0.07 

0.22 

MtlM 

3rd 

3.8*0.07 

0.22 

0.2*0.07 

2.0*0.13 

4th 

EH| 

■iW 

0.2*0.07 

5th 

m 

■iW 

wBzm 

3.8*0.08 

0.25 

Egg 

3.8*0.08 

3.8*0.08 

0.26 

0.2*0.08 

1.3*0.09 

0.09 

With  Shielding  -  No  Filter  Applied  | 

Contrast 

SNR 

Iteration 

Mean 

Noise 

Bias 

Mean 

Noise 

1th 

3.1*0.07 

0.21 

0.9*0.07 

giraiia 

0.27 

2nd 

3.5*0.07 

0.22 

0.5*0.07 

EEStEa 

0.16 

3rd 

3.6*0.07 

0.23 

0.4*0.07 

WEMM 

0.12 

BEE3B 

3.6*0.07 

0.23 

mrnnna 

1.5*0.03 

3.6*0.07 

EESI 

liSSAkl 

1.3*0.03 

■mi:* 

6th 

3.6*0.07 

ESDI 

0.4*0.07 

1.2*0.03 

7th 

3.6*0.07 

0.23 

0.4*0.07 

1.2*0.03 

8th 

3.6*0.07 

■iW 

0.4*0.07 

1.2*0.02 

The  horizontal  bars  represent  the  standard  error  of  the  mean 
calculated  by  dividing  the  noise  by  Vw  >  where  N  =10  is  the  number 
of  individual  scans. 

It  is  expected  that  with  an  OS-EM  algorithm,  the  higher 
the  iteration  number,  the  noisier  and  the  closer  to  the  true  contrast 
the  image  is.  Such  characteristics  are  shown  in  both  shielded  and 
non-shielded  cases,  although  the  effect  is  very  subtle.  To  observe 
such  behavior  more  clearly,  3-D  filtering  was  applied  to  the  20th 
iteration  of  the  reconstructions.  The  Hann  filter  as  described  below 
was  used  with  different  cutoff  (vc)  values  and  the  results  are 
tabulated  in  Table  2,  and  the  Noise  vs.  Bias  plot  is  shown  in 
Figure  4.  Since  filtering  makes  the  bias  higher,  we  expected  the 
plot  to  resemble  an  inversely  proportional  curve.  The  results  are 
mostly  consistent  with  our  predictions. 

//(v)  =  Q5»(llC0S(--)) 

VC 

From  both  Figure  3  and  Figure  4,  it  is  clear  that  shielded 
images  are  more  biased.  A  comparison  between  shielded  and  non- 
shielded  projection  views,  shown  in  Figure  5,  reveals  that  the  non- 
shielded  scans  contain  a  higher  number  of  counts  from 
approximately  180  to  260  degrees.  From  the  plot  it  is  clear  that  the 
lead  plate  is  moderately  effective  in  reducing  counts  from 
background  scattering. 


Table  2:  Shielding  and  Filtering 


Without  Shielding  -  Hann  Filter  Applied  | 

Cutoff(vc) 

Contrast 

SNR 

cycles/cm 

Mean 

Noise 

Bias 

Mean 

Noise 

1.4 

MRM1 

0.18 

■H'l 

10*0.3 

0.81 

1.6 

3.2*0.06 

0.2 

ga 

9.2±0.2 

0.72 

1.8 

3.3*0.07 

0.21 

TiTCM 

8.2*0.2 

0.63 

2 

3.4*0.07 

0.22 

»Ti 

7.3*0.2 

0.52 

2.2 

3.5*0.07 

0.22 

Bn 

6.5*0. 1 

0.44 

2.4 

3.5*0.07 

0.23 

beh 

5.8*0. 1 

0.4 

2.6 

3.6*0.07 

0.22 

5. 3*0.1 

0.36 

With  Shielding  -  Hann  Filter  Applied  | 

Contrast 

SNR 

cycles/cm 

Mean 

Noise 

Mean 

Noise 

1.4 

2.8*0.07 

0.21 

■Q 

9.2±0.3 

1 

1.6 

3*0.07 

0.22 

i 

8.4*0.3 

0.85 

1.8 

3.1*0.07 

0.24 

Baa 

7. 5*0.2 

0.71 

2 

3.2*0.08 

0.24 

EH 

6.7±0,2 

0.61 

2.2 

3.2*0.08 

0.26 

■ttEl 

6.0*0.2 

0.51 

2.4 

3.3*0.08 

0.26 

■IM 

5.4*0. 1 

0.43 

2.6 

3.3*0.08 

0.26 

beh 

4.8*0. 1 

0.37 

x  10*  Total  Counts  vs.  Projection  View 
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Figure  5:  Counts  in  Projection  Views 


However,  our  findings  in  Figure  3  and  Figure  4  are  in 
contradiction  with  those  in  Figure  5.  Reduction  in  background 
counts  (Figure  5)  is  expected  to  decrease  bias  and  increase  SNR, 
but  the  bias  and  SNR  evaluations  prove  otherwise.  The  use  of 
shielding  slightly  reduced  SNR  from  about  1.3  ±  0.09  to  1.2  ±  0.02 
on  the  f  iteration.  However,  shielding  increased  bias  in  contrast 
from  0.2  ±  0.08  to  0.4  ±  0.07.  We  believe  that  this  obvious 
discrepancy  may  be  due  to  inconsistency  in  projection  data,  or  to 
incomplete  sampling  of  the  other  organs.  As  the  gamma  camera  on 
Triad  system  rotates  around  the  phantom,  significant  difference  in 
counts  rises  when  the  camera  level  changes  from  ‘below  the  lead 
shielding’  to  ‘above  the  lead  shielding’.  With  our  current 
reconstruction  algorithm,  this  inconsistency  in  projection  data 
might  have  introduced  artifacts  that  lead  to  undesirable 
reconstruction  results. 


Conclusion: 

In  this  study,  we  investigated  the  effects  of  lead  shielding 
on  breast  imaging.  Our  results  indicate  that  shielding  is 
moderately  effective  in  reducing  the  background  counts 
suggesting  potentials  to  improve  lesion  detectibility.  However,  the 
standard  reconstruction  algorithms  such  as  OS-EM  and  ML-EM 
are  limited  in  reconstructing  inconsistent  projection  data. 
Comparison  of  contrast  and  SNR  between  shielded  and  non- 
shielded  images  revealed  that  contrast  and  SNR  is  significantly 
lower  with  lead  shielding  than  without  lead  shielding.  Future 
research  is  required  to  determine  whether  an  OS-EM  algorithm 
can  be  implemented  that  can  effectively  account  for  inconsistent 
projection  data. 
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